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We present a report on the status of the GRAL project (Going Realistic And Light), which aims at simulating 
full QCD with two dynamical Wilson quarks below the vector meson decay threshold, m ps /m v < 0.5, making use 
of finite-size-scaling techniques. 



1. INTRODUCTION 

The computer simulation of full QCD in the 
regime of light quark masses, i.e. below the vec- 
tor meson decay threshold, R = m ps /m v < 0.5, 
represents one of the most interesting and yet 
long standing goals of lattice QCD. As the algo- 
rithms and computer resources of today are still 
inadequate for directly simulating realistic quark 
masses without running into severe finite size ef- 
fects (FSE), GRAL will make a virtue of necessity 
by performing simulations on sequences of small 
and medium sized lattices; using finite size scal- 
ing techniques we will attempt to extrapolate the 
results to the infinite volume prior to the chiral 
and continuum extrapolations. In order to make 
contact with SESAM and T^L data we use stan- 
dard Wilson fermions. 

2. FIXING PHYSICAL PARAMETERS 

Seeking to simulate QCD with light quarks our 
first task is to determine an operating point (/3 op , 
k° p ) in parameter space with R° p < R d = 0.5. In 
order to render simulations at this point feasible 
we specify N s = 16 as the maximal spatial extent 
of a series of lattices employed. With the number 
of lattice points in temporal direction held fixed 
at N t = 32, the computational cost for N s = 16 
thus sets an upper bound for the costs of subse- 
quent simulations on smaller lattices. 

To be specific, we require the simulation with 
N° p = 16 at (/3 op ,K op ) to yield meson masses 
with R op = 0.4 and a FSE of z op = 0.2. 
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This particular value for the finite size parame- 
ter z = 1/ (am ps N s ) has been chosen in view of 
a recent SESAM/TxL quark mass analysis |Q, 
where z « 0.2 was found to mark the onset of 
significant FSE's in the spectrum. 

In summary, the envisaged properties of our 
target point on a 16 3 x 32-lattice are as follows: 

^ z £ ps am ps am v a[fm] L [fm] 
~QA 0.2 3.2 031 0.78 0.116 1.86 



The lattice spacing a = a p has been estimated on 
the basis of the SESAM and T%L data at /3 = 
5.6 and 5.50. Obviously, the price for sticking 
to N s — 16 is a fairly coarse lattice. However, 
according to the scaling formulae in [|j , we expect 
a total cost of approximately 0.2 TFlop/s • h to 
produce 100 statistically independent gauge field 
configurations at (/3 op , ft op ), which is feasible. 

3. SEARCHING THE TARGET POINT 

On the search for the bare run parameters 
(/3 op , k° p ) we let ourselves be guided by a heuris- 
tic consideration. From the known lattice spac- 
ings a = 0.095 fm at p = 5.5 and 0.078 fm at 5.6, 
we deduce an estimate for /3 op from the renormal- 
ization group /3-function. Assuming that we are 
in the scaling regime we get 8(3/ d In a ~ —0.46, 
which leads us to expect that the presumed lat- 
tice spacing of a = 0.116 fm at the target point 
will be attained at /3 op « 5.4. 

With this in mind we turn to the SESAM/TxL 
spectral data. The solid lines in fig. [l] represent 

1 For key parameters of the data see table 1 of 111 . 
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Figure 1. Chiral extrapolation of m v and 
for (3 = 5.6. 

linear fits to the m v and mp S data at /3 = 5.6 
with N a = 16. The dashed lines result from sep- 
arate fits to the T%L data with N s — 24. We use 
these fits (and corresponding ones for /3 = 5.5) to 
establish functional dependencies of R and z on 
1/k, as shown in fig. g. 

We are now in the position to guesstimate the 
operating point by direct extrapolation in the 
(3-k plane. Setting out from the functions of 
fig. |^ we use a linear and a quadratic ansatz 
to construct curves of constant R and z in 
that plane, displayed in fig. ||. The intersec- 
tion of the straight lines for R = 0.4 (solid) 
and z = 0.2 (dashed) in the (3-n plane delivers 
with (/3°p,k°p) = (5.32,0.1652) a first estimate 
for the operating point. Although /3 op comes out 
close to the estimate from the /3-function, this 
result is not to be taken too seriously, as pre- 
liminary results of a simulation performed with 
these action parameters give R((3 op ,n op ) = 0.8 
and z(/3 op ,K op ) = 0.1, which is off the expecta- 
tion by a factor of two in both cases. In order 
to illustrate qualitatively the behavior of i?(/3, k), 
the long-dashed curve shows a quadratic interpo- 
lation between (/3 op ,n op ) (having R = 0.8) and 
the points with the same R at (3 = 5.5 and 5.6. 
(The dotted curve shows the same for z((3, k).) 

A two-dimensional quadratic fit of all available 
(simulated and extrapolated) i?(/?, k) and z(f3, k), 
respectively, to functions of the form /(/3, n) = 
co + ci/3 + C2K, + cz/3 2 + c^k 2 + c§(3k leads to the 
coordinates (/3 op ,K op ) = (5.34,0.1653). The dot- 
dashed curves in fig. || show the resulting /3(k) 
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Figure 2. R (a) and z (b) as functions of 1/k. 

for R = 0.4 and z = 0.2. As this second guess 
for (/3 op ,K op ) deviates only little from the first 
one we are still running at (3 = 5.32, but with 
increased values of k. Production is carried out 
on APEmille at DESY Zeuthen. 

4. SIMULATION PLAN 

For the smaller lattices (N s < 16) we will use 
the cluster computer ALiCE at the University of 
Wuppertal. To carry out a chiral extrapolation, 
we will generate additional ensembles at equal /?, 
but larger R. The step to a larger j3 can be done 
either by going to N s = 18,20,22 at R = 0.4 
and z — 0.2, or by moving on to z > 0.2, depend- 
ing on the results at the present simulation point. 
We are aware of the fact that our operating point 
may be located far in the strong coupling regime. 
On the other hand we might possibly take advan- 
tage of the increased "/3-shifts" at light sea quark 
masses. Eventually we want to reach values of 
(3 > 5.6. 
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Figure 3. Lines of constant m ps /m v and z in the /3-k plane. 



For the extrapolation to the infinite volume 
limit we expect a formula of the form m = 
moo + cL~ u to be applicable S. Fig. || shows 
results of Fukugita et ai; the solid lines are fits 
with v = 3 to their full QCD data. 
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Figure 4. Extrapolation to infinite volume. 



